Statistics of charge transfer in a tunnel junction coupled to an oscillator 

J. Wabnig, D. V. KhomitskyQ and J. Rammer 
Department of Physics, Umea University, SE-901 87 Umea 
A. L. Shelankov 

Department of Physics, Umea University, SE-901 87 Umea and 

A. F. Ioffe Physico -Technical Institute, 194021 St. Petersburg, Russia. 
(Dated: 9 September 2005) 

The charge transfer statistics of a tunnel junction coupled to a quantum object is studied using 
the charge projection technique. The joint dynamics of the quantum object and the number of 
charges transferred through the junction is described by the charge specific density matrix. The 
method allows evaluating the joint probability distribution of the state of the quantum object and 
the charge state of the junction. The statistical properties of the junction current are derived from 
the charge transfer statistics using the master equation for the charge specific density matrix. The 
theory is applied to a nanoelectromechanical system, and the influence on the average current and 
the current noise of the junction is obtained for coupling to a harmonic oscillator. 
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I. INTRODUCTION 

In recent years it has become possible to couple charge 
dynamics of electrons to vibrational modes of a nanos- 
tructure, and the new field of nanoelectromechanics has 
emergedti Nanoelectromechanical devices are expected 
to lead to new technologies, such as ultra-small mass de- 
tection techniques ^ as well as stimulating fundamental 
studies of quantum phenomena in macroscopic systems^. 
For example, experiments have probed high-frequency 
nanomechanical resonators in order to reach the quan- 
tum uncertainty limit Charge transfer by mechani- 
cal motion has been studied in experiments shuttling sin- 
gle electrons iSiifliil Nanomechanical resonators are inves- 
tigated for use in quantum information processingjiSiii 

The simplest possibility of detection and control of the 
vibrational degree of freedom in a nanomechanical system 
is at present via the coupling to a quantum point contact 
or a tunnel junction. Recent theoretical studies of nano- 
electromechanical systems have therefore considered the 
coupling of a harmonic oscillator to a tunnel junction. 
The master equation for the oscillator was originally ob- 
tained by Mozyrsky and Martin^ using a method limit- 
ing considerations to the zero-temperature case, and the 
influence of the coupled oscillator on the average current 
was obtained. Smirnov, Mourokh and Horing examined 
the non-equilibrium fluctuations of an oscillator coupled 
to a biased tunnel junction^ Clerk and Girvin derived 
the master equation for the oscillator, and considered the 
current noise power spectrum in the shot noise regime, 
studying both the dc and ac cases4& Armour, Blencowe 
and Zhang considered the dynamics of a classical oscilla- 
tor coupled to a single electron transistor^! and Armour 
the induced current noised Related models were studied 
in molecular electronics. 19 



Tunnel junctions functioning as position detectors or 
being monitored by a vibrational mode put emphasis on 
a description of the current properties of a tunnel junc- 
tion in terms of its charge dynamics. Recently we con- 
sidered a general many-body system coupled to a quan- 
tum object and considered their joint dynamics in the 
charge representation^ This approach, based on a previ- 
ously introduced charge projection technique^! provides 
a quantum description of charge dynamics based directly 
on the density matrix for the system, and allows us to 
treat the number of particles in a given piece of mate- 
rial as a quantum degree of freedom, establishing thereby 
in proper quantum mechanical context the charge repre- 
sentation. The evolution of the coupled systems is de- 
scribed in terms of the charge specific density matrix for 
the quantum object, /5„(i), i.e., the dynamics conditioned 
on the number n of charges in a specified spatial region 
of the environment. When a many-body environment is 
coupled to another quantum object, the method allows 
evaluating at any moment in time the joint probability 
distribution describing the quantum state of the object 
and the number of charges in the chosen region of the 
many-body system. The charge specific density matrix 
description of the dynamics of a quantum object is there- 
fore an optimal tool to study transport in nanostructures 
since in electrical measurements any information beyond 
the charge distribution is irrelevant. So far we have ap- 
plied the method to charge counting in a tunnel junction 
coupled to a discrete quantum degree of freedom, viz. 
that of a two-level system, and shown that the charge 
state of the junction can function as a meter providing 
a projective measurement of the quantum state of the 
two-level system^ 

In this paper we shall apply the charge projection 
method to the case where the quantum object coupled to 
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the junction is a continuous degree of freedom. In partic- 
ular, we shall concentrate on the properties of the current 
through the junction due to the coupling to the quantum 
object. The statistical properties of the current through 
the junction and its correlations with the dynamics of the 
quantum object coupled to it, shall be expressed through 
the charge specific density matrix. We shall illustrate the 
results for the case of a harmonic oscillator coupled to 
the tunnel junction. It should be noticed, that current 
experimental setups studying nanoelectromechanical sys- 
tems are operated under conditions where temperature, 
oscillator excitation energy, and voltage bias across the 
junction are comparablei*^^ Furthermore, nanomechan- 
ical oscillators, such as a suspended beam, are in addition 
to the charge dynamics of the electrons in the junction 
invariably coupled to a thermal environment, say the sub- 
strate upon which the oscillator is mounted. We are thus 
considering the situation where a quantum object in ad- 
dition to interacting with a heat bath is interacting with 
an environment out of equilibrium. Having the additional 
parameter, the voltage, characterizing the environment 
in non-equilibrium, gives rise to features not present for 
an object coupled to a many-body system in equilibrium. 
The presented approach is applicable in a broad region of 
temperatures and voltages of the junction and arbitrary 
frequency of the oscillator and thus generalizes previous 
treatments. 

The paper is organized as follows: In Sec. II we in- 
troduce the model Hamiltonian for a generic electrome- 
chanical nanoresonator, a harmonic oscillator coupled to 
a tunnel junction, and derive the Markovian master equa- 
tion for the charge specific density matrix. The master 
equation for the charge unconditional density matrix, i.e., 
the charge specific density matrix traced with respect to 
the charge degree of freedom of the junction, is discussed 
in detail for the case of a harmonic oscillator coupled to 
a tunnel junction. In Sec. Ill we consider the influence 
of the oscillator on the current-voltage characteristic of 
the junction. In Sec. IV we consider the properties of the 
stationary state of the oscillator. We calculate the heat- 
ing of the oscillator due to the nonequilibrium state of 
the junction, and calculate the steady state I-V charac- 
teristic of the junction. In Sec. V we consider the current 
noise in the junction using the charge representation, and 
obtain the explicit expression for the current-current cor- 
relator in the Markovian approximation. In Sec. VI, the 
theory is then applied to the case of an oscillator influ- 
encing the current noise of the junction. Finally, in Sec. 
VII we summarize and conclude. Details of calculations 
are presented in appendices. 



II. MASTER EQUATION 

As a model of a nanoelectromechanical system we con- 
sider a harmonic oscillator coupled to a tunnel junction. 
The transparency of the tunnel barrier is assumed per- 
turbed by the displacement, x, of the oscillator. The 



resulting Hamiltonian is 

H = H + H,+H r + H T 



(2.1) 



where Hq is the Hamiltonian for the isolated harmonic 
oscillator with bare frequency and mass m. A hat 
marks operators acting on the oscillator degree of free- 
dom. The Hamiltonians Hi :T specify the isolated left and 
right electrodes of the junction 
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where l,r labels the quantum numbers of the single par- 
ticle energy eigenstates in the left and right electrodes, 
respectively, with corresponding energies e\ r . The oper- 
ator Ht describes the tunnelling, 
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with the tunneling amplitudes, Ti r = Tj,, depending on 
the oscillator degree of freedom. Due to the coupling, the 
tunnelling amplitudes and thereby the conductance of the 
tunnel junction depend on the state of the oscillator. In 
the following we assume a linear coupling between the 
oscillator position and the tunnel junction 



Wl r X 



(2.4) 



where v\ r = v*^ is the unperturbed tunneling amplitude 
and wi T = w*! its derivative with respect to the position 
of the oscillator. The derivation of the equation of mo- 
tion for the charge specific density matrix presented in 
axmendixlAl shows that the following combinations of the 
model parameters v\ r and w\ T enter the master equation: 
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(2.5) 

These lumped parameters for the junction have the fol- 
lowing physical meaning: Let G(x) = e 2 G(x), e being 
the electron charge, denote the conductance as a func- 
tion of the oscillator coordinate x when it is treated 
as a classical variable. Then, Go gives the conduc- 
tance of the junction in the absence of coupling to the 
oscillator, G = jiG\x=o, and G x = 2^rffU=o and 

G xx = ^ cb§U=o- ^ ne cou pli n g constant g x cannot be 
expressed via G(x). Note that g x changes its sign upon 
interchange of tunneling amplitudes between the states 
in the two electrodes, i.e., after the substitution 1 <-> r. 
Therefore, it is only finite for an asymmetric junction 
and is a measure of the asymmetry. As shown in section 
IIIII g x generates effects similar to charge pumping, as 
well as nontrivial features in the electric current noise as 
discussed in section IVA For later use it is convenient to 
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present the coupling constants in terms of conductances 
by introducing the characteristic length of the oscillator 



(2.6) 



where x — (h/mujo) 1 ^ 2 , and loq is the frequency of the 
coupled oscillator as introduced in appendix El 



A. Charge specific master equation 

To study the interaction of charge dynamics in a tun- 
nel junction with the dynamics of a quantum object, we 
describe the combined system, the quantum degree of 
freedom coupled to a tunnel junction, using the charge 
specific density matrix method introduced in Ref. |2ij. 
The approach employs charge projectors to study the dy- 
namics of the quantum object conditioned on the charge 
state of the junction. The charge projection operator, 
V n , projects the state of the conduction electrons in the 
junction onto its component for which exactly n electrons 
are in a given spatial region, say in the left electrode. The 
charge specific density matrix is then specified by 



p n {t) =Tr ei (V n p(t)) 



(2.7) 



where pit) is the full density matrix for the combined 
system, and Tr e ; denotes the trace with respect to the 
conduction electrons in the junction. Provided the sys- 
tem at the initial time, t — 0, is in a definite charge 
state, i.e., described by a charge specific density matrix 
of the form p n (t) = S n o Po, where po is the initial state 
of the quantum object, the charge index n can be inter- 
preted as the number of charges transferred through the 
junction. Thus the charge projector method provides a 
basis for charge counting statistics in the cases where the 
distribution function for transferred charge is relevant as 
discussed in Ref. El The charge specific density matrix 
allows therefore the evaluation, at any moment in time, 
of the joint probability of the quantum state of the ob- 
ject and the number of charges transferred through the 
junction. For example, if the charge specific density ma- 
trix is traced over the quantum object degree of freedom, 
the probability p n (t) that n charges in time span t are 
transferred through the low transparency tunnel junc- 
tion is the expectation value of the charge projector, or 
expressed in terms of the charge specific density matrix 



p„(t) = Tr(p„(i)) 



(2.8) 



where the trace is with respect to the degree of freedom 
of the coupled quantum object. 

The Markovian master equation for the charge spe- 
cific density matrix, p n {t), for the case of coupling of the 
junction to a quantum object is derived and discussed 
in appendix El The Markovian approximation is valid 
for describing slow time variations of the density ma- 
trix; the exact conditions of the applicability are spec- 
ified later once the characteristic times of the problem 



have been identified. To lowest order in the tunneling, 
the master equation for the charge specific density ma- 
trix can be generally written in terms of super-operators: 
a Lindblad-like term A, a diffusion term V and a drift 
term 



Pn 



:[H ,p n ] + A{p n } + V{&} + J{p' n }, (2.9) 



(2.10) 
(2.11) 



where p' n and p 1 ^ denote the discrete derivatives, 

Pn = \ (Pn+1 - Pn-l) , 
Pn = Pn+1 + Pn-l — %Pn- 



General expressions for the super-operators in Eq. i!2,9|) 
are presented in appendix El as well as their specific form 
for the case of coupling to an oscillator. The equation 
shall in Sec. IIVI be used to study the current noise in 
the junction due to the coupling to a quantum object, 
before we in Sec. consider the explicit case of an oscil- 
lator coupled to the junction. However, first we analyze 
the master equation for the unconditional density ma- 
trix, i.e., the charge specific density matrix traced with 
respect to the charge degree of freedom of the junction. 



B. Unconditional Master equation 

Often interest is not in the detailed information of the 
charge evolution of the tunnel junction contained in the 
charge specific density matrix. If for example interest is 
solely in properties of the oscillator, this information is 
contained in the traced charge specific density matrix. 
We are thus led to study the master equation for the 
reduced or charge unconditional density matrix, the den- 
sity matrix traced with respect to the charge degree of 
freedom, p(t) = J2 n p n (t). Performing the charge trace 
on Eq. H2.9J1 , the master equation for the reduced density 
matrix can be written in the form 

1-7 D 



A 
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(2.12) 



The form of the master equation is generic to any contin- 
uous quantum degree of freedom coupled linearly to the 
junction, and has the well-known form for a particle cou- 
pled to a heat bath i 22 i 2 i? In the following we consider the 
model Hamiltonian for a nanoelectromechanical system 
introduced in section |nj and encounter the renormalized 
oscillator Hamiltonian 



P 

2m 



(2.13) 



which in addition to having a renormalized oscillator fre- 
quency ujq, suffers a voltage dependent linear shift in the 
equilibrium position of the oscillator, which in the follow- 
ing is assumed absorbed into the position of the oscillator 
(for details see appendix EJ- 
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The second and third term on the right in Eq. 112.1211 
represent the physical influences of friction and fluctua- 
tions of the environment. For a nanomechanical object, 
the environment consists of several parts. The first one, 
which we have explicitly included in the model, is the 
tunnel junction. The other one is included phenomeno- 
logically by introducing 70 and D , the values of the fric- 
tion and diffusion parameters in the absence of coupling 
to the junction. The physical mechanism generating the 
friction coefficient 70 and diffusion coefficient D is, e.g., 
the heat exchange of the nanoscale oscillator and the bulk 
substrate it is mounted on. Thus the latter environment 
could also be modeled microscopically in the standard 
manner of coupling a quantum object to a heat bath i 22 ' 24 
Then, with the model Hamiltonian in Eq. 112. If) giving the 
electronic environment contribution due to the coupling 
to the junction, the total friction and diffusion coefficients 
become 



7 = 7o + 7e 



D = Dn 



D e 



(2.14) 



where 7 e is the electronic contribution to the damping 
coefficient 



7e 



(2.15) 



proportional to the coupling strength G xx , and the elec- 
tronic contribution to the diffusion coefficient is 



n 

D e = mj e Q coth 2^7- 



(2.16) 



where ft = huio and the voltage dependent parameter, 
T e , is given by the relation 



coth 



n 

2T P 



V 
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coth 



ft V 



2T 



2ft 



ft , V 
coth ■ 



ft 



2T 
(2.17) 

Here T is the temperature of the junction and we assume 
a dc voltage bias, V = ell, U being the applied voltage. 
We note that the r.h.s. of Eq. (I2.17|) is proportional to 
the well-known value of the power spectrum of current 
noise of the isolated junction, taken at the frequency of 
the oscillator^ The fact that the junction as a part of 
the environment is in a non-equilibrium state, is reflected 
in the voltage dependence of the electronic contribution 
to the diffusion coefficient. In section lTVl we show that T e 
is the effective temperature of the junction seen by the 
oscillator. 

The phenomenological parameters Dq and 70 are re- 
lated to each other by virtue of the fluctuation-dissipation 
theorem. Assuming that the junction and the part of the 
environment responsible for 70 and Do have the same 
temperature, T, the diffusion coefficient can be generally 
presented in the form 



D = mft 70 coth 



ft 
2T 



7 e coth 



ft 

27; 



(2.18) 



The master equation Eq. (|2.12|) contains the term pro- 
portional to the coupling constant A. A term with 



this structure has been obtained in previous discussions 
of quantum Brownian motioniS2*2£ The derivation of 
the master equation for the oscillator (see appendix El 
shows that the main contribution to the coefficient A in 
Eq. 112.1211 comes from virtual tunneling processes with 
energy difference of initial and final states of the order of 
the Fermi energy, in contrast to the friction and diffusion 
coefficients which are controlled solely by the tunneling 
events in the vicinity of the Fermi surface. Besides, com- 
pared with the other terms in Eq. Q2.12I1 . the ^4-term 
has different symmetry relative to time reversal, i.e., the 
transformation p — » (p) and t — > —t. The damping 
and diffusion terms, which are odd relative to time rever- 
sal, describe the irreversible dynamics of the oscillator, 
whereas the last term in Eq. (I2.12|) just like the Hamilto- 
nian term is time reversible. These observations give the 
hint that the A-term is responsible for renormalization- 
like effects. This suggests that the A-term should be 
treated on a different footing than the dissipative terms. 
Indirectly, the A-term can be absorbed into the Hamilto- 
nian dynamics at the price of having the time evolution 
of the oscillator described by a "renormalized" density 
matrix. Indeed, if we apply a (non-unitary) transforma- 
tion to the density matrix, p = lZ{p}, by acting on the 
density matrix p with the super-operator 



K{p} = p + p[p, [p, p], 
we can by proper choice of the parameter p, 

A 



/' 



2mft 2 



(2.19) 



(2.20) 



cancel the A-term in the equation for the transformed 
matrix p. Leaving the 7 and D— terms intact, the counter 
term is produced by the super-operator 1Z acting on the 
Hamiltonian part of the master equation Eq. <2.12t . Ap- 
plied to the original A-term, this procedure generates 
an additional contribution proportional to the product 
pA oc A 2 , and it can be neglected to lowest order in the 
coupling, the limit which can be consistently studied. 

The renormalized density matrix now obeys the master 
equation 



where 



JC{p} = -{H R ,p] 



p(t) = IC{p} 



D 



(2.21) 



x,[x,p]}, (2.22) 



up to a term quadratic in the coupling constant. For 
compact notation, we have introduced the super-operator 
/C, and dropped the tilde for marking the renormalized 
density matrix: thus in the following the renormalized 
density matrix will also be denoted by p. The master 
equation being derived for the case of coupling to a tun- 
nel junction is seen to be of the same form as for cou- 
pling linearly to a heat bath, i.e., an equilibrium state 
of a many-body system; the generic form of a damped 
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quantum oscillator known from numerous investigations 
on quantum Brownian motion^ However, the diffusion 
term is qualitatively different from the usual case where 
the quantum object is coupled only to a heat bath. The 
non-equilibrium state of the junction, characterized by 
its voltage, gives rise to features not present when the 
coupling is simply to a many-body system in equilibrium. 

We note that the super-operator, 1Z, does not change 
the trace of the density matrix it operates on, and the 
renormalized density matrix is also normalized to unity. 
However, one has to keep in mind that the observables 
are to be calculated with the "unrenormalized" density 
matrix. Up to first order in the coupling constant the 
expectation value of an observable O are now given in 
terms of the renormalized density matrix according to 

(°)=^({ d -^\P^])p). (2-23) 

This relation transfers renormalization from the density 
matrix to observables. In the language of the Feynman 
diagram technique, Eq. il.'I'M corresponds to a vertex 
correction. 

In this paper, we use the Markovian approximation to 
describe the time evolution of the density matrix. This 
approximation is valid in the low frequency limit. For the 
rfc-bias case, the characteristic frequency of time varia- 
tion of the density matrix, w, must be small enough to 
meet the condition 

(T V\ 

LO<^UJ m ax , Umax = max I — , — I . (2.24) 

From the unconditional master equation, the character- 
istic frequency is seen to be determined by the friction 
coefficient, oj ~ 7. This means that the coupling con- 
stant G xx in Eq. H2.15J1 must be small enough to meet 
the condition 7 <C Lo max ■ 

III. CURRENT- VOLTAGE CHARACTERISTIC 

The average value of the current through the junction 
is given in terms of the probability distribution for charge 
transfers, i.e.^Sl 

I(t) = -ej t J2 n ^Mt) (3.1) 

n 

where Tr denotes the trace with respect to the degree 
of freedom of the coupled quantum object. However, to 
lowest order in the tunneling, the average current turns 
out to be expressible through the charge unconditional 
density matrix, the reduced density matrix for the cou- 
pled quantum object. Indeed, the master equation for 
the charge specific density matrix then enables one to 
express the time derivative in Eq. 13. lt in terms of the re- 
duced density matrix for the quantum object, the charge 
unconditional density matrix 

I(t) =eTrJ{p{t)} (3.2) 



where the drift super-operator, J , is specified in 
Eq. {ATT)l . 

A. Contributions to the current under dc bias 

For a dc bias V = ell, U being the applied voltage, the 
drift operator J is given by Eq. l|Allll . and the current, 
Eq. H3.2H . is specified by 

~/(t) = V{G) t + I q (V) + I p (t) (3.3) 

and in general time dependent due to the coupling to 
the oscillator. The current consists of three physically 
distinct contributions. The first term is the Ohmic-like 
part of the current proportional to the conductance 

(G) t = G + 2G x (x) t +G xx (x 2 ) t , (3.4) 

the instantaneous value of the conductance operator, 
Eq. 1ET2J, where (x n ) t = Tr(x n p(t)). We note that 
besides the pure Ohmic term of the isolated junction, 
the additional terms due to the coupling to the oscillator 
will in general contribute to the non-linear part of the 
current-voltage characteristic since the state of the oscil- 
lator will depend on the voltage. A case in question is 
discussed in the next section, where the stationary state 
of the oscillator is considered. 

The second term, I q , originates according to the deriva- 
tion due to position and momentum operators are not 
commuting and for this reason we refer to it as the quan- 
tum correction to the current 

I g (V) = ~G xx A v , (3.5) 

where Ay is specified in Eq. l!A13jl or equivalently 

Ay = V + (SI + V)N n +v - (SI- V)Nq- V (3.6) 

and Nn±v = l/(e^ - 1). 
The last term in Eq. ffjfll . 

V 

I p (t) = ehg x (x) t , x = — , (3.7) 

m 

is proportional to the average velocity of the oscillator, 
and present only for an asymmetric junction, g x ^ 0. 

The Ohmic part of the current is calculated in section 
II V Bl for the stationary case. Next we turn to discuss the 
quantum correction and the dissipationless contribution 
to the current. 



B. Quantum correction to the current 

The quantum dynamics of the oscillator leads to a sup- 
presion of the Ac current as expressed by the quantum 
correction to the current, I q {V). Unlike the other terms 
in the expression for the current, Eq. i'.'i.'M . the quantum 
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FIG. 1: Voltage dependence of the quantum correction to the 
junction conductance, G q . The three curves correspond to the 
following relations between the junction temperature T and 
the oscillator frequency: T = O.Olfi, T = fi, and T — 10f2. At 
high temperatures, T fl, the quantum correction is small 
and shows no particular voltage dependence (dotted line). At 
low temperatures, T <C fi, two distinct regions can be distin- 
guished. At voltages lower than the oscillator frequency the 
quantum correction stays constant, and at voltages higher 
than the oscillator frequency it is inversely proportional to 
the voltage. 



correction, I q {V), does not depend on the state of the 
oscillator, but only on its characteristic energy and the 
temperature of the junction. At low voltages, V <C T, 
the quantum correction is linear in the voltage 



Ig = -VG X 



l+N Q --N Q (N n + l) 



(3.8) 



The presence of the term I p in the current, which does 
not depend explicitly on the voltage, means that a cur- 
rent through the junction can be induced just by the 
motion of the oscillator alone, i.e., by the time variation 
of a system parameter which in our case is the junction 
transparency. This effect is closely related to the well 
known physics of quantum pumping ^SLSi but has not, 
to our knowledge, been discussed in the present context. 
The dissipationless "pumping"- like current, I p , is propor- 
tional to the coupling constant g x . This contribution to 
the current is thus only present if the tunnel junction is 
asymmetric. This is in concordance with the quantum 
pumping effect »24 The global symmetry properties of the 
system thus crucially determine the existence and mag- 
nitude of the induced current. 

A single mode oscillator driven by an external peri- 
odic force at frequency to induces an ac-current, I~, with 
the same frequency and a phase of the ac current rigidly 
following the phase of the external force. For a given 
amplitude of the oscillations, x max , the magnitude of the 
ac-current can be estimated as 

I^~a as euj (3.10) 

where the dimensionless parameter a as — hg x x max char- 
acterizes the effective asymmetry of the junction. In prin- 
ciple, a as may be comparable to unity so that 1^ ~ eu, 
provided the amplitude of the oscillations x max is large 
enough and the conductance of the junction is not too 
small. 

One can show that in the case of an oscillator with two 
or more modes interacting with the junction, the corre- 
sponding term generates directed pumping of charged 



where Nq — l/(e^ 
constant value 

In 



1). At large voltages it reaches a 



~G M n , 7>T,fi, (3.9) 



in agreement with an earlier result obtained by a tech- 
nique valid at zero-temperature^ Our approach gener- 
alizes the expression for the current to arbitrary relations 
between junction voltage and temperature, and the fre- 
quency of the oscillator. The voltage dependence of the 
quantum correction to the conductance, G q = I q /V, is 
shown in FigQ]for different temperatures. 



C. Dissipationless current 

The last contribution in Eq. I|3.3|l to the current, I p , 
is qualitatively different from the other terms. From 
Eq. 113.711 . one sees that I p is proportional to the ve- 
locity of the oscillator. Therefore, the corresponding 
transferred charge through the junction, SQ P = I p St, is 
controlled by the coordinate of the oscillator: SQ P = 
ehg x 5(x). Being proportional to a velocity, the current 
contribution I p is odd with respect to time reversal, and 
therefore a dissipationless current. 



IV. STATIONARY STATE PROPERTIES 

In this section, we shall study the stationary state 
of the reduced density matrix for the oscillator in the 
Markovian approximation. The question arises whether 
the stationary state of the oscillator is a thermal equi- 
librium state even though the environment is in a non- 
equilibrium state as the junction is biased. According to 
Eq. 12. 2H , the stationary renormalized density matrix of 
the oscillator, p s , is determined by the equation 



IC{ Ps } = 0, 



(4.1) 



and the solution is indeed of the form of a thermal density 
matrix, p s oc exp ^—Hr/T*^ , where the temperature of 
the oscillator, T*, is specified by the relation 

n d 



coth • 



2T* 



Using Eq. H2.18fl . the temperature of the oscillator is re- 
lated to the environment temperature and the voltage 
bias according to 



coth • 



Q 

2T* 



7o 



70 + 7e 



coth • 



7e 



2T 70 + 7 e 



coth ■ 



2T e 



(4.2) 
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The average occupation number for the oscillator, N*, 
is given by the Bose function, N* — l/(e n ^ T — 1), and 
seen to be populated separately by the interaction with 
the two environments 



AT/T 



N* 



To 
7 



Nr, 



7e 
1 



(4.3) 



where N e = l/(e^ - 1)^1 

We observe, that the oscillator acquires the temper- 
ature of the bath, T* w T, if the interaction with the 
junction is weak and 70 is the dominant contribution to 
the friction, 70 ^> 7 e . The general case and the opposite 
limit where the dynamics of the junction is dominating, 
7e ^7o, we proceed to consider. 



A. Oscillator heating 

When the oscillator is well isolated and the interaction 
with the junction dominates, "f e 3> 70, the oscillator at- 
tains according to Eq. 14.21 the effective temperature of 
the junction T e , as given by Eq. II2.17J1 . As expected, in 
the absence of a bias voltage across the junction, V = 0, 
the temperature of the oscillator equals that of the junc- 
tion irrespective of its temperature. When the junction 
is biased, the oscillator is generally heated except at zero 
temperature and low voltages, and we first discuss the 
case of a junction at zero temperature. 

At zero junction temperature, T — 0, we must distin- 
guish two voltage regions. If the voltage is smaller than 
the frequency of the oscillator, V < SI, the temperature 
of the oscillator is also zero, T* — 0, independent of the 
voltage as it follows from Eq. H4.2I) . In this regime, the 
interaction with the tunneling electrons is unable to ex- 
cite the oscillator from its ground state. Heating can 
only take place beyond the voltage threshold given by 
the oscillator frequency. If instead the voltage is larger 
than the frequency of the oscillator, V > f2, the temper- 
ature of the oscillator, T* , is determined by the following 
relation to the voltage 



tanh • 



si 

2T* 



si 

W\ 



(4.4) 



At high voltages, V ^> Q, the temperature of the os- 
cillator approaches half the bias voltage, T* = V/2, in 
agreement with the result obtained in a previous study 
where the temperature of the junction was assumed to 
vanish^ 

The heating of the oscillator, its excess temperature, 
AT = T* — T, as a function of the bias, is shown in Fig.|2j 
both for the case where the coupling to the junction 
dominates and the opposite case of dominating external 
damping. The effect of the external damping is shown for 
moderate to strong external coupling, 7/7 e = 5, 10, 100. 
Increasing the coupling to the external heat bath leads 
to suppression of the heating of the oscillator, the ad- 
ditional environment acting as a heat sink. In the case 



AT/T 
4 




- 7, Ay = 0.2 
leh = 0.1 
777- = 0.01 



V/T 



FIG. 2: Heating of the oscillator, its excess temperature, 
AT = T* — T, as a function of the bias voltage for the tem- 
perature T = 0.1S7, and different ratios of 7e/7, covering from 
moderate to strong external coupling. The influence of the ex- 
ternal damping reduces the heating effect. The inset shows 
the heating with negligible external damping (7e/7 = 1) as 
a function of the bias voltage for different junction tempera- 
tures. At junction temperatures low compared to the oscilla- 
tor frequency, T <^Sl, two regions of voltage dependence can 
be distinguished with a rapid switch as the voltage passes the 
value of the frequency of the oscillator. If the bias voltage 
is smaller than the oscillator frequency, V < Q, the oscilla- 
tor temperature, T*, depends only weakly on the voltage. In 
the region where the voltage exceeds the oscillator frequency, 
V > SI, the oscillator temperature approaches V/2. 



where the coupling to the junction dominates, the inset 
shows that at low temperatures the oscillator is not ex- 
cited at voltages below the frequency of the oscillator. 

At high voltages, V ^> Cl 7 T, in the shot noise regime, 
the temperature of the oscillator T* can be found from 
Eq. H4,4jt , Just as in the case of vanishing junction tem- 
perature, the oscillator temperature approaches half the 
bias voltage, T* = V/2, at large bias, and these results 
generalize previous studies which were limited to zero 
temperature and high voltages, V ^ fl^ilfL 

In the quest for using tunnel junctions to measure the 
position of a coupled object with the ultimate precision 



set by the uncertainty principle 



5.6.7.: 



it is important to 



take into account that the measuring, involving a finite 
voltage, will invariably heat the oscillator. In this respect 
the presence of the additional heat bath, described by the 
coupling 70, is important. For example, envisioning the 
oscillator has been cooled to a temperature T much lower 
than SI and a voltage is turned on. In order to obtain an 
appreciable signal the voltage must be larger than SI. The 
oscillator will then in a time span of the order of 7 -1 be 
heated and attain a temperature for which the average 
number of quanta in the oscillator is 



N* 



7e \V\ - V 

2 7 SI 



(4.5) 



Estimating the oscillator temperature, we have T* ~ 
Sl/ln(7/7 e ). A strong environmental coupling can thus 
be beneficial for retaining the oscillator in the ground 
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state. 



B. I — V characteristic 

In the stationary state, the dc current I(V), Eq. 
is conveniently written on the form 

-I{V) = VG + \g xx ( 2VN* + (fi - V)N n -v 

- (fl + V)N n +v) ■ (4.6) 

As before, N* is the occupation number of the oscillator 
at the temperature T* . In the stationary state, only the 
coupling constant G xx is effective. 

The first term in Eq. i|4,6|) is the current through the 
isolated junction, and the remaining term describes the 
influence of the oscillator on the current. It is interesting 
to consider the latter in the limit of vanishing environ- 
ment temperature, T = 0. At zero temperature, we must 
distinguish two voltage regimes. If the voltage is smaller 
than the frequency of the oscillator, V < £1, we observed 
in the previous section that the oscillator remains in the 
ground state, N* — 0, and the average current turns out 
to be equal to that of an isolated junction, / = eGoV. 
One observes that the effect of zero point fluctuations 
present in the conductance, Eq. (13.41 . is exactly can- 
celled by the quantum correction, Eq. H3.5J) . This lack 
of influence of zero point fluctuations is expected since 
the oscillator in its ground state is inert to the tunneling 
electrons for such low bias. 

If the voltage is larger than the oscillator frequency, 
V > ft, the slope of the I-V characteristic, at T = 0, 
abruptly increases, / = VGq + (V — Cl)G xx /2, 

At arbitrary temperatures, the linear conductance of 
the junction, G = I/V, V — » 0, is given by 

±-G = G Q + G xx % N n (N n + l) . (4.7) 

To derive this formula, we recall that in the limit of van- 
ishing bias, V — > 0, the oscillator attains the temperature 
of the junction, T* = T. 

At high temperatures, T » V, ft, the quantum correc- 
tion to the current, and thereby the nonlinear quantum 
corrections, vanish, and we obtain the result 

-I(V) = V(G + G xx (x 2 )*) , {x 2 )* = -^ (4.8) 

where (x 2 )* is the mean square of the oscillator coordi- 
nate at temperature T*. This result is to be expected 
from a classical oscillator in thermal equilibrium influ- 
encing the conductance of a tunnel junction. The I —V 
characteristic is in this regime nonlinear due to the volt- 
age dependence of T* , Eq. 14.21 . 



V. CURRENT NOISE 

In this section we use the charge projection technique 
to develop the description of the statistical properties of 
the current of a tunnel junction coupled to a quantum ob- 
ject. The discussion will be kept quite general before we 
in the next section specialize to the case of a harmonic 
oscillator, the nano-electromechanical model of section 
Irfl We shall show how the charge dynamics, described 
by the master equation for the charge specific density 
matrix, can be used to obtain the statistical properties 
of the junction current, such as the noise power spec- 
trum. The prerequisite for the success of this endeavour 
is that for the considered low transparency tunnel junc- 
tion, the charge representation in fact provides the prob- 
ability distribution for the charges transferred through 
the junction al 21 



A. Current noise in the charge representation 

The probability, p n {t), for n charge-transfers in time 
span t is according to Eq. (I2.8|l given by p n (t) = 
Tip n (t), where p n (t) is the charge specific density ma- 
trix, Eq. The charge-transfer probability distribu- 
tion specifies the stochastic process of charge transfers, 
n(t). The variance of the charge fluctuations 

((n 2 (t))) = (n 2 (t))-(n(t)) 2 , (5.1) 

is defined in terms of the moments of the probability 
distribution of charge transfers 

(n r (t)) = ]T n r p n (t) , Pn (t) = Tr (£,(*)) . (5.2) 

n 

To express the statistical properties of the current in 
terms of the probabilities for charge transfers, we inherit 
the stochastic current process, i(t), through its relation 
to the charge transfer process 

n(t) = [ dt'i(t'). (5.3) 
Jo 

The average current, given by (i(t)) — d(n(t)) /dt, is 
in accordance with Eq. Ij.S.lJl ■ The variance of the charge 
fluctuations are expressed via the current fluctuations ac- 
cording to 

<(™ 2 W)H / dt i I #2 (Si{h)6i(t 2 )) . (5.4) 
Jo Jo 

where 5i(t) = i(t) — (i), and (i) is the average dc current 
since we in the following shall consider the stationary 
state. Stationary current noise is characterized by the 
current-current correlator 

S(t) = (5i(t + r)5i{t)) . (5.5) 
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Inserting this expression into Eq. llo.4ll one obtains 

((n\t)))=2 f dr{t-T)S{r). (5.6) 



This expression allows one to relate the charge and cur- 
rent fluctuations. 

Taking time derivatives of Eq. Il5,6|) gives 



d{{n\t))) 



and 



dt 



S(t) 



2 / drS(r), 



ld 2 «n 2 (i)» 



dt 2 



(5.7) 



(5.8) 



i.e., the current-current correlator equals the second 
derivative of the variance of charge transfers. This re- 
lation allows one to calculate the current-current correla- 
tor, S(t), by evaluating the charge fluctuations using the 
master equation for the charge specific density matrix. 

Eventually, interest is in the current noise power spec- 
trum, Su, given by 



S u = 4 / dt cos(tut)S(t) 
Jo 



(5.9) 



where u is the frequency at which the noise is measured^ 
We observe that the zero frequency noise power, ac- 
cording to Eqs. 115.711 and 115. 9J1 . can be calculated from 
the general relation 



Su>=0 



d((n 2 (t))) 



dt 



(5.10) 



i.e., as the rate of change of the charge variance at large 
times. 

In the present approach it is convenient to calculate di- 
rectly the current-current correlator, as done in the next 
section and appendix El Only at the end we then trans- 
form to obtain the noise power spectrum. However, we 
note that the approach is equivalent to employing the 
widely used MacDonald formula^ 



B. Current-current correlator 

In this section we show how the master equation for 
the charge specific density matrix, can be used to obtain 
the noise power spectrum of the current. A convenient 
feature of the method is that it allows one directly to 
obtain the time dependence of the current noise. 

The probability distribution of charge transfers, p n (t), 
is obtained from the master equation for the charge 
specific density matrix given the initial condition corre- 
sponding to a state of definite initial charge 



at the time when the charge counting starts, the initial 
time t = 0. We are interested in the noise properties of 
the stationary state and, therefore, the stationary density 
matrix of the oscillator, the thermal state p s , enters the 
initial condition. 

In the following we shall treat the charge specific dy- 
namics in the Markovian approximation. The charge 
specific density matrix p n (t), is obtained as the solution 
of the master equation Eq. I|2.9H . For notational conve- 
nience we write the charge specific master equation in 
the form 



dpn y. r * , 



V{p'n} + J{Pn}- 



(5.12) 



where /C is the super-operator introduced in Eq. 12. 2H . 
Although the p„'s are time dependent, the unconditional 
density matrix, p = ^2 n Pn{t), remains equal to the 
thermal state, p s , by virtue of its stationarity property, 
K-{f>s\ — 0. Our goal is now to evaluate the variance 
of the charge transfers and thereby the current-current 
correlator with the help of Eq. 115. 8JI . 

The rate of change of the first charge moment, i.e., 
the dc current according to Eq. l.'t.H . becomes in the 
stationary state 



1 



I = -Hi (J {p.}), 



(5.13) 



following from Eq. Ij3.2l) and the stationarity property 
of the density matrix for the coupled quantum object, 
p(t) = p s . The dc current was calculated in section HV Bl 
It readily follows Eq. H5.12J1 . that the time derivative 
of the variance of charge transfers, ((n 2 (i))), can be pre- 
sented in the form 



j t {{n 2 (t))) = 2Tr (V {p s }) - 2Tr (j {dN(t)}) 



where 6N(t) denotes the traceless matrix 



(5.14) 



= X> -<»(*)» *.(*)• ( 5 - 15 ) 

71 

We observe that only the truncated density matrix, 
SN(t), is needed to calculate the noise. 

Comparing Eq. H5.14J1 and Eq. H5.7I) . one concludes 
that the current-current correlator has a S— function like 
singularity at the initial time, t — 0, where the charge 
counting starts. Indeed, the r.h.s. of Eq. H5.14J1 has a 
finite limit as t — ► 0, given by the first term, since the 
second term initially vanishes, SN(t = 0) = 0. For this 
result to be compatible with Eq. I|5.7H and Eq. I|5.8ll . the 
current-current correlator, S(t), must have the following 
structure 



S(t) = Si(t) + S 2 (t) 
the sum of a singular contribution, 



(5.16) 



Pn(t = 0) = S nfi p a 



(5.11) 



S 1 (t)=2Tr(V{p s }) S(t) 



(5.17) 
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where S(t) denotes a function peaked at t = and nor- 
malized according to the condition f Q dt S(t) = |, and 
a regular part given by 



VI. 



NOISE POWER SPECTRUM 



S 2 (t) 



-Tr 



j{si(t)} 



(5.18) 



where SI denotes the matrix, SI = 4? SN. The finite time 
correlation of the current described by the regular part 
52 (t) is solely due to the interaction with the quantum 
object, as follows from SI(t) being traceless. We note 
here that the <5-function singularity, which would provide 
noise at arbitrary high frequencies, is an artefact of the 
Markovian approximation. 

The task of calculating the time dependent current 
noise is thus reduced to obtaining the time derivative 
of the charge-averaged density matrix, SN(t), given in 
Eq. I|5.15|l . From the master equation for the charge spe- 
cific density matrix one obtains the following equation 
for <JI(t), 



dt 



SI = JC{8I}, 



and the initial condition 



51 



t=o 



-SJ{p s }. 



(5.20) 



Here the super-operator SJ acts on its argument matrix 
according to 



SJ {X} = J {X} - X (Tr J {X}) . 



(5.21) 



We note, that acting on a matrix X with unit trace, 
TrX = 1, the super-operator SJ returns a traceless ma- 
trix. The dynamics of the charge averaged quantity SI is 
thus identical to that of the charge unconditional density 
matrix of the oscillator. 

The formal solution to Eq. H5.19J1 can be written in 
terms of the time evolution super-operator for the charge 
unconditional density matrix of the oscillator 



Hi 



(5.22) 



SI{t) = -U t {5J{p s }} 



(5.23) 

and the regular part of the current-current correlator can 
be written on the form 



S 2 {t) = Tr (SJ {U t {SJ {p s }}}) . (5.24) 

Here J in Eq. 115,1 8J1 has been replaced for SJ in 
Eq. H5.18J1 : the replacement is valid under the trace op- 
eration since Tr 5I(t) = 0. Combined with the singular 
part in Eq. Q5.16I1 . this gives the general expression in the 
Markovian approximation for the current-current corre- 
lator of a tunnel junction interacting with a quantum 
system in its stationary state. The current noise correla- 
tor has thus conveniently been written with the help of 
the Markovian super-operators IC, J and V. 

In the next section we shall turn to calculating the 
noise properties for the case of the nanoelectromechanical 
device described in section ITfl 



We now turn to calculate the current-current correlator 
of the tunnel junction coupled to the harmonic oscillator 
as described by the model of section^ Taking advantage 
of the general analysis in the Markovian approximation 
presented above in Section FV Bl the current-current cor- 
relator can be written in the form: 



S(t) = S!(t) + S x (t) + S x ,(t) . 



(6.1) 



Here Si is the singular part defined in Eq. i!5,17JI and 
specified in Eq. ijBllll . The regular contribution is given 
by second and third terms on the right in Eq. H6.ll) . as 
obtained by inserting the expression for the drift super- 
operator J, Eq. (IA11I . into the expression for the regular 
contribution, Eq. H5.24I) . giving 



S x (t) = 2VG X xj(t) + g x — pj(t) (6.2) 

m 



(5.19) and 



S x 2(t) = VG xx x 2 j(t) (6.3) 
where the time dependent quantities are given by 



Xj(t) =Tr [XUt{SJ{p s }}) , 



(6.4) 



with X = x,p,x 2 , respectively. It is readily checked that 
the latter quantities evolve in time in accordance with 
their corresponding classical equations of motion for a 
damped oscillator with initial conditions given by 



Xj{Q) = Tr(XSJ{p s } 



X = x,p, x 2 ,p 2 . (6.5) 



The calculation of these quantities are presented in ap- 
pendix El giving according to Eq. id. 21 and Eq. I|6.3ll an 
explicit expression for the current-current correlator S(t) 
and thereby, according to Eq. H5.9I) . for the noise power 
Su. Fourier transforming the current-current correlator 
gives, according to Eqs. l|B91 IB 101 IB 1 111 . peaks in the 
noise power as well as a constant up-shift in the noise 
floor, the noise pedestal. The noise power spectrum is 
displayed in the inset in Fig. 03 As expected, there is a 
pronounced peak at the frequency of the oscillator, and 
in addition two side peaks each shifted by the frequency 
of the oscillator, one at zero frequency and one at twice 
the oscillator frequency. 

Below, we analyze the noise in two frequency regions: 
(i) low frequency noise at frequencies uu < 7; and (ii) 
noise in the vicinity of the oscillator resonance frequency 
lo ps ujq and lo ps 2lu<j. We examine the voltage and 
temperature features of the noise power. 



A. Low frequency noise 

Let us consider low frequency noise, at frequencies of 
the order of the damping rate and lower, u> < 7. Then the 
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noise power spectrum is given by the Fourier transform of 
the correlation functions Si(t) and S x 2(t) in Eqs. HBlfljl . 
and HBlljl . respectively, giving 

5. = + <?« + (6.6) 

where 

5(°)=2G ^coth^ (6.7) 

is the low frequency, w <C V, white Nyquist or Schottky 
noise of the isolated junction, and is the correction to 
the white noise due to the interaction with the oscillator 

= 2G xx Sl (N* (N e + 1) + N e (N* + 1)) . (6.8) 

Together these two contributions form the noise pedestal, 
Soo = + S^. The frequency dependent part, 
becomes at low frequencies 

S { J } = G XX J™ 2 ^ (2VN* 2 + (V — A v ) (2N* + 1)) . 

(6.9) 

The low frequency noise is displayed explicitly propor- 
tional to the coupling to the electronic tunnel junction 
environment as we have taken advantage of the relation 
7e = G xx Sl. The width of the low frequency peak is 
twice the damping rate 2j. At zero bias, V = 0, where 

(2) 

Su = and the oscillator and effective junction tem- 
peratures equal the environment temperature, leaving 
N* = N e — Nq, one recovers the fluctuation-dissipation 
relation for the noise power, Sw=o = 4TG, where G is the 
linear conductance of the junction in the presence of the 
interaction with the oscillator, i.e., given by Eq. 114.711 . 

In the following we discuss the features of the low fre- 
quency excess noise, the noise due to the coupling to the 
oscillator, and in particular the noise peak height at zero 
frequency, in the limits of temperatures high and low 
compared to the oscillator frequency. 



At voltages much higher than threshold, \V\ 3> SI, the 
peak height relative to the pedestal becomes 

The zero frequency noise is proportional to V 4 if the ef- 
fective coupling to the electronic environment is appre- 
ciable, i.e., the ratio 7 e /7 is not too small. In the high 
voltage limit, the oscillator is in the classical regime but 
with the oscillator temperature given by T* = \V\/2 as 
discussed in section ITvl 

2. High temperature noise 

At temperatures higher than the oscillator frequency, 
T » SI, we can distinguish two voltage regimes. At low 
voltages, 7 <T, the peak height scales quadratically in 
both the temperature and voltage 

S%l Q = ~Gl x V*(^\ (6.12) 

and we recall that the oscillator temperature equals the 
junction temperature, T* = T. 

At high voltages, V 3> T, the peak height becomes 

S u=0 --V G xx \ -- + -—J , (6.13) 

At high temperatures, the oscillator is in the classi- 
cal regime and the average occupation number depends 
linearly on the oscillator temperature. The peak height, 
proportional to the fluctuations in the oscillator position 
squared, is proportional to the square of the average oc- 
cupation number, and is therefore proportional to the 
square of the oscillator temperature. 

B. High frequency noise 



1. Low temperature noise 

First we consider the low frequency noise at low tem- 
peratures, T <C fl. As expected, no excess noise is ac- 
cording to Eq. i|6.6Jl generated by the oscillator at zero 
temperature and voltages below the oscillator frequency, 
|V| < SI, where the oscillator cannot be excited from its 
ground state. Indeed, in the region of low temperatures, 
T < O, and low voltages, V < ft, the oscillator is non- 
responsive and the excess noise, and \ vanishes 
exponentially in O/T below the activation energy SI. 

Close to the noise onset threshold, |V| ~ SI, in the 
narrow region, ||V| — Sl\ «T, the peak height relative to 
the pedestal rises linearly with temperature 

S ( Jio = ~TG XX . (6.10) 



Next, we investigate the properties of the peaks in 
the noise power spectrum occurring at finite frequen- 
cies, at the oscillator frequency, uo « ljq, and its har- 
monic, uj an 2ujq. The Markovian approximation allows 
us to consider the high frequency noise only under the 
condition Eq. 12.241 . that the frequency is much smaller 
than the maximum value of the voltage or the temper- 
ature, and for frequencies in question, this means that 
max(T, V) ^> SI for consistency. The inset in Figure 
shows the frequency dependence of the noise power spec- 
trum, Eq. 15.9|l . in the case of high temperatures, T ^ SI. 
The noise power displays three peaks. The noise power 
spectrum, consists at uj pa ujq of a Lorentzian part, as 
given by the Fourier transform of Eq. l|B9t , with a width 
given by the damping rate, 7, and an asymmetric part 
specified by Eq. IIB9J1 and present only for an asymmet- 
ric junction. At oj w 2loq, the noise power spectrum is 
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according to Eq. (IB 10p a Lorentzian with a width given 
by twice the damping rate, 27. We now turn to discuss 
these peak heights at high and low temperatures. 



1. High frequency noise at high temperatures 

At high temperatures, T 3> V 3> O, the height of the 
peak at the oscillator frequency relative to the pedestal 
depends linearly on temperature 



S (2) - 2- 



~ 2 4y 2 _ 2 



(6.14) 



and is determined by the conductances G x and g x . 

At double the oscillator frequency the peak height de- 
pends quadratically on temperature 



C(2) 



2u.'(i 



—&c, 

7 



V 2 T 2 

~n 2 



(6.15) 



and just as the peak at zero frequency determined by the 
conductance G xx . We note that its height is half that of 
the peak at zero frequency, Eq. Ijfi.1211 . 

The expressions for the excess noise, Eq. I|6.14|l and 
Eq. Q6.15I1 . are in fact also valid at low voltage. Contrary 
to the peaks at w = and u = 2uj , which vanish in the 
absence of voltage, the excess noise power at u ~ uj is 
therefore finite for an asymmetric junction even at zero 
voltage, and according to Eq. I|6.14ll in fact negative. An 
asymmetric junction with a Q-factor much larger than 
T/O can thus at zero voltage lead to a suppression of the 
noise power below that of an isolated junction. 



2. High frequency noise at low temperatures 

At high voltages and low temperatures, V 3> fl » T, 
the peak height at the oscillator frequency becomes 



c(2) 



2GiV 1 



jfl 1- 



7e 



Je\V[ 

7 J 3X " \~ 2 7 , 

(6.16) 

and the peak height at twice the oscillator frequency 



2V 

7 



o(2) ^1^2 V 2 7e\V\_f 1 e\Vl 

UJ=2uJ0 7 xx 7 20 U 20 



(6.17) 



The noise can be large due to the high oscillator temper- 
ature. 



3. Noise asymmetry 

A striking feature of the finite frequency noise is the 
contribution proportional to g x G x . It is odd relative to 
the sign of the voltage and does not depend on the state 
of the oscillator (see Eq. IIBlflJl ). This term, which is 
only present for an asymmetric junction, g x 7^ 0, does 
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FIG. 3: The excess noise power spectrum, AS^ = — Sao, 
due to an oscillator coupled to an asymmetric junction is 
shown by the full line for the parameter values V/Q = 10, 
T/fi = 0.01, 70/wo — 0.05 and conductances Go = 1, 
G x = 0.07, g x = 0.05, and G xx = 0.005 giving for the elec- 
tronic coupling constant 7e/wo = 0.005. The dotted line 
shows the contribution symmetric in the voltage, ,and 
the dash-dotted line the contribution asymmetric in the volt- 
age, ■ The inset shows the noise power spectrum for the 
high temperature T = 10057, but otherwise the same set of 
parameters, a regime where the peaks at zero and twice the 
oscillator frequency are also visible. 



not contribute to the peak height at the oscillator fre- 
quency, but provides the asymmetry of the peak in the 
frequency region around the oscillator frequency, u> = luq. 
Separating the even and odd voltage contributions in the 
noise power, = \ (S U (V) ± S u (— V)), the odd contri- 
bution becomes 



S- = g x G x (Fuj (w) - -F_a, (w)) 
I V 2 coth — + QV(2N e + 



(6.18) 



A 

2 Av 



with the frequency dependence given by the function 

2(w - oj) 



7 2 + (lj - ui Q ) 2 



The noise power spectrum is displayed by the full line 
in Figure for the temperature T = 0.01O, where only 
the peak at the oscillator frequency is appreciable. The 
even part in voltage of the noise power, 5+, displayed by 
the dotted line, is a symmetric function of the frequency 
relative to the frequency of the oscillator, and the odd 
part in the voltage, S~ , displayed by the dash-dotted line, 
is an anti-symmetric function of the frequency relative to 
the oscillator frequency. If the voltage is reversed, the 
frequency dependence of the asymmetric part is mirrored 
around the frequency lo = ujq. 

In contrast to an isolated junction, the noise power of 
the coupled junction-oscillator system shows asymmetry 
in the voltage. This behavior is a novel feature that arises 
when an asymmetric junction is coupled to an additional 
degree of freedom. 
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VII. CONCLUSIONS 

We have applied the charge projection technique to ob- 
tain the charge specific dynamics of a continuous quan- 
tum degree of freedom coupled to a tunnel junction. The 
master equation for the charge specific density matrix has 
been derived, describing the charge conditioned dynam- 
ics of the coupled object as well as the charge transfer- 
statistics of the junction. The method allows evaluating 
at any moment in time the joint probability distribution 
describing the quantum state of the object and the num- 
ber of charges transferred through the junction. 

The approach, generally valid for any quantum object 
coupled to the junction, has been applied to the generic 
case of a nanoelectromechanical system, a harmonic oscil- 
lator coupled to the charge dynamics of a tunnel junction. 
In this regard it is important that the method allows in- 
clusion of a thermal environment in addition to the elec- 
tronic environment of the tunnel junction since nanores- 
onators are invariably coupled to a substrate. The os- 
cillator dynamics, described by the reduced density ma- 
trix for the harmonic oscillator, the charge specific den- 
sity matrix traced with respect to the charge index, has 
upon a renormalization been shown to satisfy a mas- 
ter equation of the generic form valid for coupling to 
a heat bath. Even though the electronic environment 
is in a non-equilibrium state, the master equation is of 
the Caldeira-Leggett type, consisting of a damping and 
a fluctuation term. Though the coefficients of the terms 
are not related by the equilibrium fluctuation-dissipation 
relation, the fluctuation term originating from the cou- 
pling to the junction is of the steady-state fluctuation- 
dissipation type, containing the current noise power spec- 
tral function of the isolated junction taken at the fre- 
quency of the oscillator. The diffusion parameter is thus 
determined by all energy scales of the problem including 
temperature, voltage and oscillator frequency. The pres- 
ence of an environment in a non-equilibrium state thus 
leads to features which are absent when the oscillator is 
only coupled to a heat bath. 

The Markovian master equation for the charge specific 
density matrix has been used to calculate the current. 
In general, the average junction current consists of an 
Ohmic term, however, with a conductance modified due 
to the coupling to the oscillator dynamics, a quantum 
correction, and a dissipationless ac current only present 
for an asymmetric junction and proportional to the in- 
stantaneous velocity of the oscillator. The latter term 
does not depend on voltage explicitly and is an example 
of an effect similar to quantum pumping. 

The stationary state of the oscillator has been shown 
to be a thermal state even though the environment is 
in a non-equilibrium state. Therefore, the only effect of 
the bias is heating of the junction. Thus the stationary 
oscillator state is a thermal equilibrium state, though in 
equilibrium at a higher temperature than that of the en- 
vironment if the junction is in a non-equilibrium state of 
finite voltage. This is a backaction effect of the measur- 



ing device, the tunnel junction, on the oscillator. At zero 
temperature and voltages below the oscillator frequency, 
the oscillator remains, to lowest order in the tunneling, 
in its ground state, and the dc current equals that of 
an isolated junction. The coupling of the oscillator to 
the additional heat bath, described by the coupling con- 
stant 70, is shown to be beneficial for avoiding heating 
of the oscillator due to a finite voltage. This is of im- 
portance for application of quantum point contacts and 
tunnel junctions to position measurements aiming at a 
precision reaching the quantum limit. 

The charge projection method has been used to infer 
the statistical properties of the junction current from the 
charge probability distribution. For example, the noise 
power spectrum is specified in terms of the variance of the 
charge distribution. The master equation for the charge 
specific density matrix can therefore be used to obtain 
the current-current correlator directly, and this has been 
done explicitly in the Markovian approximation. The 
excess noise power spectrum due to the coupling to the 
oscillator consists of a main peak located at the oscilla- 
tor frequency and two smaller peaks located at zero fre- 
quency and twice the oscillator frequency, respectively. 
The peaks at zero frequency and at twice the oscillator 
frequency have heights proportional to the coupling con- 
stant G xx squared, whereas the height of the peak at the 
oscillator frequency is proportional to the coupling con- 
stants G x and g x squared. The voltage and temperature 
dependencies of the peaks has been examined in detail. 

For an asymmetric junction, the noise power spectrum 
contains a term with the striking feature of being an odd 
function of the voltage and independent of the state of the 
oscillator. Contrary to the case of a symmetric junction, 
coupling of an oscillator to an asymmetric junction with 
temperature higher than the oscillator frequency results 
even at zero voltage in a suppression of the noise power 
at the oscillator frequency, the excess noise power being 
negative. For an asymmetric junction, the noise power at 
uj « ljq can thus be suppressed below the Nyquist level 
of the isolated junction. 

The Markovian approximation employed to calculate 
the noise power can not be validated at arbitrary fre- 
quencies compared to temperature or voltage. Not sur- 
prisingly, naive attempts to extend expressions beyond 
the Markovian applicability range, Eq. 112. 2411 . leads to 
unphysical results for the noise. For example, at zero 
temperature and voltages below the oscillator frequency, 
a spurious noise power arises even for the oscillator in the 
ground state. 
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APPENDIX A: CHARGE SPECIFIC MASTER 
EQUATION 

In a previous paper we introduced the charge represen- 
tation for a general many-body system^ The approach 
is based on the use of charge projectors previously in- 
troduced in the context of counting statistics i21 In the 
charge representation, the dynamics of a quantum ob- 
ject coupled to a many-body system is described by the 
charge specific density matrix 



p n {t) = Tr e l (V„ p(t)) 



(Al) 



where p(t) is the full density matrix for a many-body 
system and a quantum object coupled to it, and Tr e ; de- 
notes the trace with respect to the degrees of freedom 
of the many-body system, in the following assumed the 
conduction electrons of a tunnel junction. The charge 
projection operators V n , which project the state of the 
system onto its component for which exactly n electrons 
are in a specified region of space, have been discussed in 
detail earlier i 20 i 21 There we discussed the circumstances 
under which the charge index, n, can be interpreted as 
the number of charges transferred through the junction, 
and the charge projector method thus provides a basis for 
charge counting statistics in the cases where the distribu- 
tion function for transferred charge is a relevant concept. 
In this case, the charge specific density matrix allows the 
evaluation, at any moment in time, of the joint proba- 
bility of the quantum state of the object and the num- 
ber of charges transferred through the junction. In the 
previous paper, the non-Markovian master equation for 
the charge specific density matrix for an arbitrary quan- 
tum object coupled to a low transparency tunnel junction 
was derived^ A non-Markovian master equation is less 
tractable for calculational purposes and the Markovian 
approximation is employed in the present paper. This 
is quite sufficient for calculations of average properties, 
such as the average current through the tunnel junction, 
where only the long time behavior needs to be addressed. 
However, when calculating the current noise of the junc- 
tion, the Markovian approximation limits the description 
to the low frequency noise as discussed in section IvTl 

The Markovian charge specific master equation for a 
quantum object coupled to the junction was in general 
shown to have the form^fi 

Kit) = hH ,p n (t)]+A{p n (t)}+v{p';(t)}+j{pUt)}, 

in 

(A2) 

where p' n and p" denote the "discrete derivatives" intro- 
duced in Eqs. 112.101 l2~TTl) . and the Lindblad-like super 
operator, A{p n }, has the form 



aw = s 



£>(i-A)( 

lr 



T, 



ft 



5) (A3) 



r 



1 \r 



Ti, 



H.c. 



where here and in the following H.c. represents the her- 
mitian conjugate term with respect to the variable of the 
quantum object. The bracket denotes the operation 



dT ^(e+eV)T/h e -iH T/hf ire tH T/h^ 



where T\ r is the oscillator perturbed tunneling amplitude, 
Eq. H2.4H . and e = 6\ — e r , and f\ and / r are the single 
particle energy distribution functions for the electrodes 
which in the following are assumed in equilibrium de- 
scribed by the junction temperature T. In this paper, 
we restrict ourself to the case where the junction is bi- 
ased by a constant voltage U, denoting V — ell, where 
e is the electron charge. The dagger indicates hermitian 
conjugation of operators of the coupled quantum object. 
The drift super-operator is 



J{R} = 



lr 

f ([« 



r, 



RT\ r + T\ r R 



RT\ r — T\ r R 



T, 



TV 



H.c 



H.c 



(A5) 



and it has been written in terms of the symmetric and 
antisymmetric combinations of the distribution functions 

K = h + St - 2/l/r, K = h - It. 

The diffusion super-operator can be obtained from the 
drift super-operator according to 



1 



T3{R} - - J s „ a {R}, 



(A6) 



where the subscript indicates that symmetric and an- 
tisymmetric combinations of the distribution functions 
should be interchanged, Ff T <-> F^ r . 

The bracket notation is specified in Eq. l|A4ll . where 
in general Hq denotes the Hamiltonian for the isolated 
arbitrary quantum object. In the following we consider 
an oscillator coupled to the junction, and Hq represents 
the isolated harmonic oscillator. 

As expected, the coupling of the oscillator to the tun- 
nel junction leads to a renormalization of its frequency, 
ujg — > u>q. The renormalization originates technically 
in the term present in the Lindblad-like operator, A{p}, 
which is quadratic in the oscillator coordinate and of 
commutator form with the charge specific density ma- 
trix, and gives for the renormalized frequency 

3» 1 ° r (P + + P_)|u, lr | 2 



lr 



where 



f*4 



e - V =F ftwo + *0 



(A7) 



For the considered interaction, the renormalization can 
be simply handled by changing in Eq. l|A4ll from the evo- 
lution by the bare oscillator Hamiltonian to the oscillator 
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Hamiltonian with the renormalized frequency, the shift Hamiltonian, Hq, in 1IA4JI the renormalized Hamiltonian, 

being compensated by subtracting an identical counter Eq. 112. 131) . all quantities are expressed in terms of the 

term. The above frequency shift is then identified by physically observed oscillator frequency ojo. In particu- 

the counter term having to cancel the quadratic oscil- lar, the bracket becomes 
lator term of commutator form generated by the A{/5}- 
part in Eq. (IA2j) . Substituting for the bare oscillator 









-A 



2 



2mwn 



Wlr 

2 



r 



where 



6 =S(e-V), S± = S(e-VTtiio ), (A9) 



In the following the notation — hu>o for the character- 
istic oscillator energy is introduced. 



and 



1 



e-V + iO 



Evaluating the diffusion and drift operators for the 
case of position coupling of the oscillator to the junction, 
(A10) Eq. H2.4H . we obtain for the drift super-operator 



J{R} = Vg{R} + ^G xx A v ^(xRp) + G x ^A v [p,R]+g x ^{p,R} + g x ^ (vcoth^ + Sy ) >-. />'. (All) 



r 



where the conductance super-operator is defined as the latter being proportional to the current noise power 

spectrum at the frequency of the oscillator. For the dif- 
0{R} = GqR + G x {x, R} + G xx xRx, (A12) fusion super-operator we obtain 



and 

v + n , v + n v-ci , v-n 

Ay = — - — coth — — z — coth ■ 



2T 



2T 



(A13) 



and 

v + n , v + n v-n , v-n 

Sv = —j- coth + —j- coth 



(A14) 



V{R} = - coth ^Q{R} + ^f- (b v xRx + ^(xRp)j + ^R(pRx) 



(A15) 



r 



where 



and we have introduced the notation 



B v = S v - V coth ^ (A16) K(pRx ) = i ({pRx) + (pRx)^ , (A17) 
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and 



^(pRx 



) = 1 (fpi&) - (pRx?) , (A18) 



in ijAllj) and (IA15fl . 

The parameter ^4 in Eq. IIA15JI is given by 



-4, 



2mVL 



J2\mA 2 F? r (P+-P-), 



(A19) 



and was encountered and discussed in connection with 
the unconditional master equation, Eq. 12.121 Techni- 
cally, it originates in our model from the principal value 
of integrals, i.e., from virtual processes where electronic 
states far from the Fermi surface are involved. Estimat- 
ing its magnitude under the assumption that the cou- 
plings |u>i r | 2 are constants, one obtains with logarithmic 
accuracy 



2f\? G t.t . 
A w — In 



E F 



max(V,T,0) 



(A20) 



In the course of evaluating the diffusion and drift op- 
erators, combinations like {v{ r w\ r ) appear together 
with principal value terms. The phase of vf r w\ r will in 
general be a random function of the electron reservoir 
quantum numbers 1 and r. Summing over these quan- 
tum numbers, where the principal value term does not 
provide any restriction of the energy interval, as it hap- 
pens in the case of terms proportional to delta functions, 
they will tend to average to zero, and we shall therefore 
in the following neglect such termsi^ 



APPENDIX B: NOISE IN MARKOVIAN 
APPROXIMATION 

In this section we evaluate in the Markovian approxi- 
mation the current-current correlator of the tunnel junc- 
tion for the case of a harmonic oscillator coupled to the 
junction. The task has been reduced to evaluating the 
expressions in Eq. I|6.2|) and Eq. 116. i.e., quantities of 
the form Eq. 116. 4J1 where the involved super-operator is 
the evolution operator for the charge unconditional den- 
sity matrix given in Eq. ll2.21ll . 

It immediately follows from the master equation, 
Eq. il2~2l")l . that quantities like X(t) = Tr(Xp(t)), where 
X can denote x, p, x 2 , and p 2 , and p is an arbitrarily 



normalized solution to the master equation, satisfy the 
corresponding classical equations of motion for a damped 
oscillator. The variables entering Eq. 116. 4J1 can therefore 
be expressed in terms of their initial values at time t = 0, 
and restricting ourselves for simplicity to the case of weak 
damping, 7 <C uq, they have the form corresponding to 
that of an underdamped classical oscillator 



Xj(t) = xj(0)e 7 *cosw i + 



and 

Pj(t) 
and 



Pj(0) . 



e" 7t sinw £ (Bl) 



-mu>oXj(0)e ^ s'mu> t + pj(0)e 7 *coswoi 

(B2) 



Xj(0)cOS UlQt 



p 2 (0) 



,2,.,?' 



sin toot 



(B3) 



The initial values in these equations are found from 
Eq. lE3 to be 



(B4) 



Pj(0) 



= g x - Ul(2N* + l)- \ Vcoth— + S v 



(B5) 



P 2 A0) 
2m 



G xx I SIVN +Cl(V — Ay) (iV* + — ) ) (B6) 



mu 2 x 2 (0) p 2 (0) 



2m 



G xx QVN*(N* + 1) (B7) 



{x,p}j(0) =0 



(B8) 



Substituting these initial values into equations Eqs. (IBlt - 
IB3J1 . we obtain using Eqs. I|6.2|l and (|6.3I) the expressions 
in Eq. Ij6.ll) for the regular part of the current-current 
correlator 



S x (t) = Gle'^ cos uj tV {2V(2N* + 1) - Ay) - ~g 2 x e^ 1 cos uj t QftFcoth ^ + Q 2 (N e - N*)^j 

-g x G x e'^ sincjoi (v 2 coth ^ + nV{2N e + 1) - V K j (B9) 
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and 

S x 2(t) = -Gl x e~ 2nlt V (2VN* 2 + (V- A v ) (2N* + 1)) + G^e'^V 2 N* (N* + l)cos2w t . (BIO) 

I 



Evaluating in Eq. Q5.17JI the trace of the diffusion super- 
operator in the stationary state of the oscillator, V{p s }, 



We observe that as to be expected, the Markovian 
approximation for the dynamics of the charge specific 
density matrix only captures the low-frequency noise, 
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into the integrand in the expression for the power spec- 
trum, Eq. 15.9(1 . Upon inserting the expression for the 
current-current correlator, Eq. 15.811 . into the expression for 
the power spectrum, a partial integration gives MacDon- 

ald's formula 5 W = 2ui J Q °° dt e~ 6t sin(ut) d ^" dt ^^ , since 
the vanishing of the boundary term at t = is ensured by 
virtue of Eq. 15.71 . 

A similar renormalization was encountered and discussed 
in connection with a two-level system coupled to a tunnel 



junction— 

We recall the result for an isolated junction! 20 using the 
non-Markovian master equation for the charge specific den- 
sity matrix, one obtains for the power spectrum of current 
noise = I(V + to) coth + I(V — lo) coth the 
general result for a biased junction, valid for arbitrary fre- 
quency and non-linear I — V characteristic, I — I(V). The 
Markovian approximation only captures the low frequency 
limit of the noise power. 



